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ABSTRACT 

Presented paper describes the method for finding the intersection of class space 
rational Bezier curves. The problem curve/curve intersection belongs among basic 
geometric problems and the aim of this article is to describe the new technique to 
solve the problem using relief perspective and Bezier clipping. 

Keywords: Bezier curve, perspective collineation, Bezier clipping 



1 Introduction 

Curve/curve intersection is one of the fundamental problems of computational 
geometry. At the present time there exist several different approaches to this prob- 
lem but the endeavor is to avoid difficulties in calculations which are mainly results 
of polynomial representation rational higher degree curves. It can be done in cer- 
tain cases by using the relief perspective and Bezier clipping. 



Bezier clipping, in the context of plane curve intersection in this paper, see e.g. 
Nishi90|| |[Nishi98|| , is an interactive method which takes advantages of the convex 



hull property of Bezier curves. Regions of one curve which are guaranteed do not 
intersect a second curve can be identified and subdivided away. 



Relief perspective is a mapping of space into space in order to correspond 
conditions of the human seeing. Images of all objects under the relief perspective 
are located in space between two parallel planes. 



2 Relief perspective 

The relief perspective is a special case of perspective collineation of the ex- 
tended Euclidean space E 3 . We remind that a collineation is a bijective mapping 
tp : E 3 — > E 3 that preserves collinearity of points. If there exists a plane u such 
that f{X) = X MX G u, if is a perspective collineation (for more details, see e.g. 
||Bus53|| ||Cizm84|| ). Then there exists a point O (called the centre of perspective 



collineation) such that <p(0) = O and p(a) = a Wa; O E a. 



The relief perspective is a perspective collineation including some additional 
conditions: 

• in order to produce correct images of 3D objects, it is important to respect 
necessary conditions of the human seeing. As in linear perspective, see e.g. 



|Cenek59| |, also in the relief perspective we suppose, that objects are located 
inside viewing circular cone. The cone has a vertex in the eye (the centre of 
projection) and the distance from the eye to the objects has to be at least 
25 cm. 

• the image plane 10 (the set of invariant points ) does not contain the centre 
O (i.e. <p is a homology) and determine elements of the mapping p are: the 
centre O, the image plane a (the set of invariant points) and the vanishing 
plane uf (the image of the plane at infinity). 

• no objects are ideal, i.e. objects and their images in the relief perspective 
have not ideal points. In notions of E 3 the image plane has to be placed 
between the point O and the plane u r . The mapped object, that image we 
want to construct, is located in the semi-space, that it is opposite to the 

semi-space aO ("behind the image plane a"). 

Let p be a relief perspective of space E 3 . Let denote E 3 be a preimage space and 
E3 be an image space (E 3 with upper index "r") such that ip : E 3 — > Eg. The image 
of the object WeE 3 under the relief perspective <p is called the relief of the object 
U and is denoted analogically W '. The relief perspective is given by the centre O, 
the image plane a and ordered pair of different points A, A r = <p(A) (A, <p(A) 7^ O) 
such that A, A r a and O, A, A r are collinear. 

Note therefore all mapped objects U and their images U r under the relief per- 
spective have no ideal points, we can study them in E 3 or E 3 using affine and 
projective methods as well. It is also clear that the restriction p/U of the relief 
perspective (p, which is a bijective mapping, is again a bijection. 



3 Representation of the relief perspective in an analytical 
form 

Suppose X and X' are two different points of E 3 C E 3 with coordinates [x, y, z] 
and [x', y', z'\ respectively. If p is a collineation of the space E 3 we use the following 
equations: 

a\\X + a 12 y + a>i 3 z + a M 



/ 

x 



a^x + a^y + + 



, a 2 \x + a 22 y + a 23 z + a 24 M , 

y = (i) 

a^x + a^ 2 y + a^z + 044 

, a 3i x + a 32 y + a 33 z + a 34 

z = , 

a 4 ix + a^ 2 y + a 43 z + a 44 
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where ajj G {1..4}) are real numbers and det(a,ij) 7^ 0. In equations ([!]) are 
used non-homogenic coordinates, which correspond to a restriction (p = y>/{E 3 — 
a}, where a is a plane 041X1 + 042X2 + O43X3 + a 44 xo = 0. 

In order to obtain simple mapping equations for relief perspective from ([!]) we 
assume, that the centre O is a point with coordinates [0, 0, 0] and a is a plane 
z — 1 — 0. In this special case we obtain the convenient representation of the relief 
perspective in the following form 



x' 



[1 + k)x 
z + k 



z' 



[l + k)z 
z + k 



where the span k G M + — {1} (the geometric representation of the parameter k 
will be explained in the next lines). 



4 Some properties of the relief perspective 

• the vanishing plane u r = 1 + k and the neutral plane v = —k (the preimage 
of the plane at infinity) are parallel and for distances O, v and uj t , a we have 

I Ov | = | uj r a |= k 

where the parameter k represents the span 

• if a plane a = z — c (c 7^ 0) is parallel to the image plane (a 7^ a, v), then 
the relief of the plane a is the plane a r = and these planes are parallel 
(it is obvious that a \\ a r \\ a) 

• the relief of a point placed in the semi-space, that is opposite to the semi- 
space <rO, is a point in a part of space with boundary planes z = 1 and 
z = 1 + k (in an intersection of spaces z > 1 and z < 1 + k) (Fig. |l|) 

• the orthographic projection of the object U r (the relief of the object U) onto 
the image plane a is the central projection of the object U from the point S, 
where S £ I n v (l;0 € I A I _L a) 



5 Rational Bezier curves and relief perspective 

The relief perspective is a mapping of the extended Euclidean space E 3 . Let 
Vi, i G {0,1, ... ,n} be given points, [xj, y$] be their coordinates of the space E 2 and 
positive real numbers Wi, i G {0, 1, . . . , n} be their weights. By these elements and 
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Figure 1 



Bernstein polynomials B?(t), see e.g. ||Farin93| , is defined a n-th planar rational 
Bezier curve 



n n 



P(*) = m ], *e<o,i> (3) 

i=0 i=0 

and also the space nonrational Bezier curve 

n n n 

P(t) = \^x iWi B?{t), ^ViWiBW), J2^B?(t)], t G (0, 1) (4) 

i=0 i=0 i=0 

with control points 

Vi = [xiWt, yiWi, Wi]. (5) 

Let ip : E 3 — > Eg be the relief perspective of the space E 3 and X r [x, y, z] (X r £ a) 
be the relief of the point X[a, b, c] G E 3 . From the equations @ we get 

(l + k)a (l + k)b (l + k)c 

X = ; , V = ; — , Z = : 6 

c + k ' y c + k c + k v ; 

and according to the relation z = ^jr^ we obtain 

. k(l + k) 
c + k 



1 + k-z 



From the above results coordinates [a, b, c] of the point X, that is the preimage of 
the relief X r [x, y, z], are computed according to a = x ^_ z , b = 1+ ^_ z , c = 1+ ^f_ z ■ 



Lemma Let (p : E 3 — > E 3 5e relief perspective of the space E 3 . T/ze point [x, y, z] 

fcx ky kz 

-1+k-z > 1+k-z > 1+k- 



is the relief of the point r — ^ — 
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By given planar points Vi and their weights Wi, where 1 + k > W{ Vz, we can 
define in E 3 points 

^ = t , J > i i !, ; T~p z£ {0,l,...,n} 7 

1 + K — Wj 1 + K — Wj 1 + k — Wi 

With respect to the lemma above it is obvious, that the reliefs of the points Wi are 
represented by the points (||). Using the points V*, the space nonrational Bezier 
curve P(t) of the form (§) is defined and we notice, that P(t) under mapping ip is 
an image of the curve 

n n n 

E kx iWi B^(t) £ ky iWi B?{t) E fcw f B?(i) 
Q(*) = br^ . ir^ , ] (8) 

E(i + E(i + Eii + t-^t) 

i=0 i=0 i=0 

This curve Q(t) is the space rational Bezier curve defined by control points of the 
form (|7|) and their weights Qi, that are computed according to fij = 1 + /c — Wi 
because of 

kxiWi kyiWi kwi 



1 + k — Wi 1 + k — Wi ' 1 + k — Wi 

Applying all these results we obtain the following theorem 

Theorem The curve ~P(t) is a relief of the curve Q(t). 



In order to get more information about Bezier curves and the relief per- 
spective, let us assume having a central projection with the centre O[0, 0,0] and 
the plane z — 1 = 0. In this case the planar rational curve P(t) of the form (||) is 
the image of the space rational Bezier curve Q(t) defined by (||). 

Let the space nonrational Bezier curve P(£) expressed by @ be given. What 

n 

is the relief of this curve? We know z + k = EC^i + k)-^f CO an d now the relief 

i=0 

P r (t) of the given curve P r (t) can be written as 

n n n 

J2(i + k)x iWi B?(t) E(i + %i^?(0 E(i + 

P F (0 = i—n , —n , ^ ] (9) 

EK + WCO EK + fc)£?(0 EK + k ) B Ut) 

i=0 i=0 i=0 

£ G (0, 1), where control vertices are points expressed as follow 

y. T _ (1 + /Qa;^ (1 + k)yjWi (1 + k)wj 
Wi + k ' + k ' Wi + k 

and weights f2, are computed by = Wi + k. 
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6 Intersection of space rational Bezier curves 

Let P(t), Q(u) be space rational Bezier curves of the form (§). The aim is to 
find their intersection. 

The solution to this problem can be found using already known facts. Ac- 
cording to the theorem, which was formed in the previous section, reliefs of space 
rational Bezier curves of the form (||) are space nonrational Bezier curves of the 
form (HD (in Fig. Bezier curves of the form (||) and their reliefs for k = 2 are 
shown). The central projections of these nonrational curves from the point [0, 0, 0] 




Figure 2 

onto the plane z — 1 = are planar rational Bezier curves defined by (|). Their 
intersection is possible to find specifically by using the method of Bezier clipping. 
This method specifies values of the parameter t G (0, 1) of the curve P s (t), respec- 
tively the parameter u G (0, 1) of the curve Q s (w), that correspond to the common 
point R s (the upper indexes "r" and "s" denote the reliefs and central projections 
of curves or points). This point is the intersection of both curves . The preimage 
of the point R s in the central projection is the space point, which is the relief of 
the intersection of the space rational Bezier curves P(t) and Q(u). 

The following scheme shows a whole proces of finding the intersection of the 
P(t) and Q(u) curves: 

p( t ) pr (t ) JU p-(t) 1 B .^ mg Rg 

Q(u) Q r («) Q s («) J 

R s = p s (t) n Q s (u) ^R r ^R 

R = P(t) n Q(u) 

if - relief perspektive : E 3 — > E 3 
ip - central projection ip : E 3 — > E 2 . 
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The central projection ip is not a bijective mapping and in case that the inter- 
section of curves P s (t) and Q s (w) exists (in opposite case the given curves P(£) 
and Q(w) certainly do not intersect) the intersection of curves P(t) and Q(w) does 
not have to exist. This situation occurs when preimages of point R s on the curves 
P r (t) and Q r (w) in the central projection are different. 

7 Conclusions and future work 

The relations between Bezier curves and the relief perspective have been de- 
scribed. The necessary and sufficient conditions for expressing the space nonra- 
tional Bezier curve as the relief of the space rational Bezier curve have been formed. 
We have shown that to the planar rational curve of the form @, defined by planar 
points Vi G E 2 and their weights tOj, it is possible to assign a class of the space 
curves by the relief perspective. One of them is nonrational curve defined by (||) 
and two are rational curves defined by (19) and (^. 

The described method can be used as a direction for application how to find the 
intersection of the space rational curves of the form @. It is possible to express 
every polynomial curve in Bezier's representation and due to this representation 
the method can be applied to all polynomial curves after modifications (e.g. spline 
curves which are considered as curves consist of Bezier segments) for solution to 
the curve/curve or curve/line intersection problems. 

Despite the author's attempt he did not succeed in finding similar comparable 
published methods on website. In addition to this he is not able to compare his 
results with any other. Obviously author's knowledge is limited and he would 
appreciate to get information about any other similar method. 

In the future work we want to extend possibilities of the relief perspective in 
geometric modeling. Our aim is to find an answer whether any polynomial 3D 
curve can be converted to the space curve of the form @. 
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